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PREFACE 


. Emphasis continues to ba placed on the usa of digital computers in 
solving nonlinear hydrodynamic and thermodynamic equations of fluid 
flow. This publication of the thermal pollution group at the University 
of Miami presents the solution of one such problem. This problem deals 
with the use of a numerical one-dimensional model in predicting the tem- 
perature profiles of a deep body of water. Although this model can be 
applied to most lakes, a specific site (Lake Keowee, S. C.) application 
has been chosen and described in detail. The programs are written in 
fortran V and could be modified by the user. Some of theso modification# 
are suggested either in the text or in the specific programs. 

* 

A detailed derivation of the equations integrated has been left out; 
however, to improve readability of the final equations, the meaning of 
the terms and variables occurring in these equations are included. 

This research was performed at the thermal pollution laboratory at 
the University of Miami. Funding was provided by the National Aeronau- 
tics and Space Administration (NASA-KSC) and the Environmental Pro- 
tection Agency (EPA-RTP) . 
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ABSTRACT 


A user's manual for a one-dimensional thermal model Is described. 
The model is essentially a set of partial differential equations which are 
solved by finite difference methods using a high speed digital computer. 
The main equations integrated are discussed. The programs are written 
In fortran V and an example problem Is discussed in detail. 
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SECTION 1 
INTRODUCTION 


It Is important that tha tharmal behavior of heatad dischargas and 
thair receiving basins be clearly understood. 

A numerical model that can be used for predicting tha seasonal thermo- 
clina of a deep body of water is vary useful in studying the environmental 
impact of thermal dischargas from power plants. This is not only required 
for existing power plants but also for planned units. Thus, a predictive 
capability Is essential to the licensing procedure. Monitoring programs 
cannot sa'tlsfy these needs, but from time to time, play a vital role in the 
calibration and verification of mathematical models. 

The one-dimensional, thermal numerical model, described in this 
manual, features the effects of area change with depth, nonlinear inter- 
action of wind-generated turbulence and buoyancy, absorption of radiative 
heat flux below the surface, thermal discharges and the effects of vertical 
convection caused by discharge. The main assumption in the formulation 
of this model is horizontal homogeneity. 

This model can be applied to most stratified deep bodies of water. 

This stratification has a seasonal cycle and is an important natural 
characteristic of a body of water. The body of water could b? divided 
into any number of slices. The temperature of each slice is predicted 
by the model. The surface slice exchanges heat with the environment 
of known climatic conditions while the bottom slice is assumed perfectly 
insulated. Condenser cooling water is extracted from any one of the 
slices and heated by the power plant. The discharge is injected into a 
slice of the same temperature as the dischage. 

The main function of the model is the prediction of the temperature 
profiles in a deep body of water for any number of annual cycles. How- 
ever, predictions connot be made on hourly basis - a feature usually 
handled by a more sensitive three-dimensional model. This is tha main 
limitation of the model. 

The procedure used in writing this manual is as follows: 

Description and flow chart of the main program are given in Section 
3, where the subroutines are also described. In the next section, a 
list of the variables and dimensions are given. The next three sections 


t 


* 

show how a typical run is prapartd, axacutad and plotted. An example 
case Is discussed in Appendix A, while Appendix B gives the fortran 
source program listings. 
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SECTION 2 


RECOMMENDATIONS 


The main disadvantage of a one-dimensional thermal modal Has In 
tha fact that resolution is sacrificed for computational spaed. Three di- 
mensional models are bulky and time consuming but have much better 
resolution, however, when long term simulations are necessary, a one- 
dimensional model Is recommended. 

The model described here can be modified to Include the single 
effects of the various quantities involved in the surface heat transfer 
phenomenon rather than using the equilibrium temperature concept. 

This is particularly recommended for the user who is interested In 
modeling the long term effects of one (for example, evaporation) of 
the quantities involved in the surface heat transfer processes.* 

Furthermore, Jthe modal can be easily adapted to handle connected 
multiple domains. This recommendation Is discussed in- the text. 
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SECTION 3 

PROGRAM DESCRIPTION AND FLOW CHART 


DESCRIPTION OF PROGRAM ALGORITHM 
Background 

A view of an Idealized deep body of water is shown In Figure 1* 

This basin Is divided Into eleven slices. The inner nine slices are of 
«*qual thickness. DZ, while the top and bottom slices are of thickness 
DZ/2. The thickness, DZ, Is determined from the depth of the basin 
and the number of slices used. The temperature of each slice is as 
shown in Figure 1;' the horizontal lines correspond to the center of each 
slice. 

The condenser cooling water (CCW), if any, could be taken from 
any slice. In Figure 1, the CCW is extracted from the center of Slice 2 
which is at temperature T-. The discharge temperature, T Q , is the sum 
of T, and the Increase in temperature through the condenser. T Q is In- 
jected into a slice of equal temperature or treated as a surface outfall if 
T D is greater than the highest temperature of the basin. 

The basin also gains or loses heat from the surface as a result of 
changing climatic conditions which are required as input data. These 
could vary every time step, daily or monthly. 

Algorithm 

The problem Is an initial value problem, so the values of dependent 
variables are assumed known initially. The governing and associated 
equations are discussed in the next section. The governing equation is 
parabolic and mathematically represents a diffusion process with vertical 
convection. 

The values of the dependent variables at successive time steps are 
obtained by using a forward-time Dufort-Frankei scheme. 

The sequence in which calculations are performed is as follows: 
(Refer to Summary of Variables - next section.) 

1. The dependent variables, T, K?, W*, Aw, p, T~ and Ke* are ini- 
tialized. The area of each slice is calculated ana then tne time step 
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Figure 1. Idealized deep body of water 


5 




\$ calculated. . Tha heading of tha beginning year Is printed. The 
values of the variables, K-, W*, A v , p, T £ and K-, are then calcu- 
lated. The temperatures 6T the slices are finally calculated, if the 
temperature profile Is unstable, mixing of the unstable portion of 
the profile is undertaken. 

2. During the. next time step, the temperatures are updated, and the 
dependent variables are calculated again. 

3. The values of the temperature T, eddy diffusivity K z , number of 
days and surface heat transfer coefficient K- are prfhted every time 
step, every day or normally at the end of each month. At the end 
of the present year, the title of the new year is printed and compu- 
tions continue as listed above. These steps are shown In a flow chart, 
Figure 2. The results are stored on a magnetic tape and plotted 
when necessary. 

Description of Main and Subprograms 

* 

The fortran calculation programs consist of a main program (NASA) 

and seven subroutines (YEARS, EQUIL1, STORE, CCW, SMOOTH, MIXIT 

and AREAS). 

1. MAIN; The main program handles the input data, calls the subrou- 
tines and does the temperature calculations. Two alternatives are 
given for handling the input data; these are either read through 
cards or in-data files or through a block-data arrangement given 

at the beginning of the main program. For users interested in the 
block-data package, the following caution is necessary: Whenever 
a data or set of data is changed, the main program must be recom- 
piled! 

2. YEARS: This subroutine prints the year heading. It is called at 
the beginning of a new year. 

3. EQUILI? This subroutine reads the dewpoint temperature, wind 
speed and solar radiation. It then computes the surface heat trans- 
fer coefficient and the equilibrium temperature. Depending on how 
the data has been averaged (e.g. days, months or years); it is 
called as often as needed. 

4. STORE: This subroutine stores the calculated data on magnetic 
tape designated as Unit 8. The stored data could be read by the 
plotting subroutine called READER. This subroutine and other plot 
programs are described later. 

5. CCW: This subroutine supplies the condenser cooling water data. 

The data is also converted to the required units by this subroutine. 

6. SMOOTH: This subroutine finds the largest value of the eddy dif- 
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Figure 2. Flow chart (calculation) 
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fusivity and uses it to calculate tha variable time step, it also 
smoothens the calculated eddy diffvislvlty for unstable temperature 

gradients. It is called every time step. 

♦ 

7. MIXIT: This subroutine looks for unstable temperature gradients 
and mixes or stabilizes the temperatures, it is also called every 
time rtep. 

8. AREAS: This subroutine handles the surface areas of each slice 
and converts the values to the required units. It is called only 
once at the beginning of the computations. 

9. INPUT : This is an in~data element containing all input data. 
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SECTION 4 

DESCRIPTION OF PROGRAM SYMBOLS 


Introduction 


The programs have been written to Calculate, as a function of depth, 
thermal diffusivity and temperature profiles over complete annual cycles. 
The equation integrated is 

A (Z)ft (pC p T) =ft fpc P A tzi K zli> - fz fpc p A «) TV z ) + ° A ' + A (Z) H (Z> 

( 1 ) 

The above equation requires two boundary conditions and one ititial con- 
dition. 


The initial condition is an input quantity supplied by the user and 
equals the homothermai temperature of the basin. The boundary condi- 
tions are: 


1. At the surface; 


k- aT 

Z3Z |Z=h 


= K s (T e 



( 2 ) 


Where Z = vertical coordinate measured from the deepest point 
Tp = equilibrium temperature 
t 5 = surface temperature 
K| a surface heat exchange coefficient 

2. At the bottom; 

Perfect insulation is assumed. 


3T 

7T |Z=0 


0 


(3) 


Calculations of the temperature profiles are made by numerical integration 
of Equation (1). Calculations start with the homothermai conditions and 
a forward explicit scheme is used. 


Each time step, the surface temperature, Tg = T12, is calculated 
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and then the temperature of each slice is calculated. Solar radiation is 
absorbed at the surface slice and the unabsorbed portion is transmitted 
exponentially to the slices below. 

The empirical relations involved in this manual are summarized below. 
A full discussion is given in the final report, Lee et al. (1980). 

Description of Main Variables 

1. Density, p, fortran variable - ROW: 

p = A 1 + B t T + C^T 1 


(4) 


where A. = density at 0°C 
1 » 1.02943 gm/cc 
B , = constant 
= - 0.00002 
= constant 
= -0.0000048 


’1 


’1 


2. Eddy diffusivity, K^, fortran variable = XKZ 

K z = ^20^ + °1*V”'* 

and 


(5) 


R, = 


_ V’z* 


aT 


( 6 ) 


i - IT 

where R. = Richardson number 

crJ = 0.1, an empirical constant, fortran variable - SIGMA 
g 1 = acceleration due to gravity, fortran variable - G 
W*= friction velocity, fortran variable - FRVEL 

= <Vp) 


= Aj + Bj(T - <t) + Cj(T - K)> (6a) 

fortran variable for o^,. A V 

where A- = 0, volumetric coefficient of expansion at 4°C, fortran 
variable - Al 


B- =* constant, fqrtran variable - A 2 
= 1.538 x 10“* 

C 7 - constant, fortran variable - A3 
= -2.037 x 10“' 


ay can also be estimated by using Equation (4) 
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where K,., = eddy diffuslvity under neutral condition (varies with 
20 time), fortran variable - XKZO 

k Zq * A 3 + B 3 sln(||gt + C 3 ) (t is in days) (6b) 

where 


3. Heat source, H, fortran variable - F6 

H * n (1 - B)Aj Z j * 0 exp(-n(Z - h)) (7) 

where 0 * 0. 5, fraction of the solar radiation absorbed at the surface 
n m 0.75, solar radiation absorption coefficient 

= net solar radiation reaching the water surface (input 
„ 0 variable) , fortran variable - HSOL 


A, * average value of K 7n , fortran 
B 3 » half annual varlatloirof K ZQ ,' f 
Cj » phase angle, foreran variame - 


fortran variable - R9 
K--,‘ fortran variable 
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SECTION 5 


a 

PREPARATION OF INPUT DATA 


The input data is storad in an in-data file - INPUT. Alternatlvaly, 
it could ba punched on cards. Tha input data is read in with an opan 
format. The main variables read are: dewpoint temperature, wind speed 

and solar radiation. In some cases where the dewpoint temperature is 
not available, the relative humidity, air temperature and a pschometric 
chart are used to find the dewpoint temperature, if this involves a lot 
of chart reading, subroutine EQUIL1 could be modified and the dewpoint 
temperature calculated from a known equation supplied by the user. 

If the latter case is used, then the input data base is enlarged to read 
air temperature, relative humidity, wind speed and solar radiation. A 
detailed input list of the constants is given in Appendix A. • 
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SECTION 6 


PLOTTINC PROGRAMS AND EXECUTION ELEMENTS 


DESCRIPTION OF PROGRAMS 

The fortran plotting routine consists of one main program (PLOTTER) 
and one subroutine (READER). 

PLOTTER: This program calls the calcomp fortran subroutines (re- 

fer to a Calcomp plotting manual for details) and the subroutine (READER) 
which reads the calculated results from a magnetic tape designated as 
Unit 8. (See Item A. 4.) A flow chart is shown in Figure 3. 

READER: Reads the calculated data stored on Unit 8 (magnetic 
tape) . 

Execution Elements 

Two execution elements are used, one for executing the calculated 
results and the other for executing the plots. 

OO-IT : This element compiles and prints the main program (NASA) 
and then prepares an entry point table, maps the necessary programs 
and subprograms, calls the In-data element containing the input data 
and finally, executes the calculations. This is done as follows for a 
UNIVAC 1100 computer at the University of Miami. 

Only one magnetic tape Is necessary. 

1. 0 ASG, AX FILE. 

The 'FILE' is assigned for the run. 

2. 0 ASC, T 8., 16N, TAPENAME 

A magnetic tape file named '8.' is being assigned. The tape is 9- 
track, and the reel number is 'TAPENAME'. The calculated results 
are stored on this tape. 

t 

3. 9 PRT, S FILS. NASA 

The main program is printed. 


13 



4. 9 PACK PILE. 

The 'FILE 1 is packed, 

5. 9 PREP FILE. 

The entry point table is prepared, 

6. 9 MAP, S 

7. IN FILE. NASA 

8. LIB FILE. 

9. END 

10. 9 XQT 

11. 9 ADD FILE. INPUT 

12. 9 FIN 

PLOT-IT: Similar to DO-IT, bu^ handles the plotting executions. 

For a UNIVAC 1100 computer the following cards are necessary. Two 
magnetic tapes are necessary. 

1. 9 ASG, AX FILE. 

2. 9 ASG, T 8., 16N, TAPENAME 

3. 9 ASG, T 11., 16, PLOTTAPE 

* 

A magnetic tape file named 'll.' is being assigned. The tape is 7- 
track , and the reel number is 'PLOTTAPE 1 . The plots are stored on 
this tape. 

4. 9 PRT, S FILE. PLOTTER 
The plot program is printed. 

5. 9 PACK FILE. 

6. 9 PREP FILE. 

7. 9 MAP, S 

8. IN FILE. PLOTTER 

9. LIB FILE. 

* 

10, END 
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11. 9 XQT 

12. 9 ADD FILE. INPUT 

13. 9 FIN 
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Figure 3. Flow chart (plots) 

* 
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APPENDIX A 
EXAMPLE PROBLEM 


The modal described iri this manual was variflad using monthly-* 
avaragtd data supplied by Duke Power Company for Lake Keowee, South 
Carolina. Accordingly, the data discussed below apply to Lake Keowee. 

SITE DESCRIPTION 

Lake Keowee is located 40 km west of Greenville, South Carolina. 

It Is the source of cooling water for Oconee Nuclear Station (ONS). It 
was formed from 1968 through 1971 by damming the Little and Keowee 
rivers. A connecting canal (maximum depth 30. 5 m) joins the two main 
arms of the lake. Flow out of the take is through the Keowee Hydro 
Station. Lake Keowee also exchanges water with Lake Jocassce-pumped 
storage station. The three-unit ONS with a net capacity of 2580 Mwe 
started operating in July 1973. ONS operated on annual gross thermal 
capacity factors of 11, 28, 69 and 59% in the years 1973 through 1976, 
respectively. From 1977 to 1979 the factors varied from 65 to 75%. A 
map showing the geometry of the lake is given in Figure 4. 

PROBLEM STATEMENT 

Calculation of Parameters and Input Data 

1. The fortran variable DM (I, J) is a two-dimensional array containing 
the temperatures at the connecting channel between Lake Keowee and 
the Jocassee-pumped storage station. The data is averaged monthly. 
The units are in degrees Celclu- (°C). I is the year counter and J 
is the month counter. The inp.-s for the first year are punched on 
the first card, the next year on the next card, and so on. Accord- 
ingly, each card contains twelve inputs in open format (real floating 
point numbers). 

2. The following fortran variables /constants are also read in with open 
format, five on one card. 

ZYEAR : starting year - 1 971 (could be changed) . 

t 

DZ: thickness of an inner slice (ft) - (maximum depth of lake) / ( 1 0. 0) . 

XKZL: lower limit of the eddy diffusivity (ft J /day) - corresponds to 
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Figure 4. Lake Keowee 




the thermal diffusivlty of solid water (IS ft 3 /day). 

H: maximum depth of lake, ft (150 ft). 

C: acceleration due to gravity (ft/sec 3 ). 

PI: ir = 3.1413926. 

A1: corresponds to A2 in Equation (6a); Al = 0 °C -1 . 

A2: corresponds to B2 in Equation (6a); A2 = 1.538 x 10~ 5 °C. 

A3: corresponds to C2 in Equation (6a); A3 = -2.037 x 10* 7 °C. 

A4: corresponds to At in Equation (4); A4 = 1.02943 gm/cc °C. 

AS: corresponds to B1 in Equation (4); A5 = 0.00002 gm/cc °C. 

A6: .corresponds to Cl in Equation (4); A6 m -0.0000048 gm/cc °C 3 . 

(NOTE: The units for A4 through A6 are automatically converted to 
consistent units in the main program.) 

TO: homothermal temperature of lake (initial condition); TO = 7.8 °C. 

C,.: specific heat; C =1.8 BTU/'lb °C. 

P P 

SIGMA: see Equation (5); SIGMA = =0.1. 

**R6,R7,R8: the friction velocities ( r / p ) are calculated for the whole 

period and fitted into a s sine curve: (friction velocity 

OMEGA) 

W* = R6 + R7 sin(|~time + R8) 

where R6 = average value of W*. 0.1 ft /sec. 

R7 = average value of the half annual variations 
of W*. 0.0 2S ft /sec. 

R8 = phase angle. 2.61 radians 
TIME is in days, not specified. 

R8.R9.R10: correspond to C3, A3, and B 3 of Equation (6b) respec- 

tively; R9 a 800 ft 3 /day and R10 = 200 ft 3 /day. 

DATA1: 0 or 1 (see below). 

3. The next set of inputs is the dewpoint temperatures, wind. speed and 


^Alternatively, friction velocity could be read in as monthly averages. 
If this alternative is followed, then DATA! = 1, otherwise DATA! = 0. 
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solar radiation. Thasa can aithar ba punched on cards or stored in 
an in-data element. They are read every month. Each card contains 
three members. For example: for January-March 1971 (Lake Keowee) , 

the data are 


3.0, 6.69, 167.0 


0., 9.3, 264.4 
6. 3, 9. 28, 264.4 

The first number on each line (each card) is the dewpoint tempera- 
ture in °C. The second one is the wind speed in ft/sec. The third 
quantity is the solar radiation in BTU/ft J day. If DATA1 = 1, a fourth 
number must be included on each line (every card) . This fourth 
quantity is the computed friction velocity for each month. 

NOTE: The in-data element described above is called INPUT,. (See 

* Fortran Source Program Listing, Appendix B.) 
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Sample Output and Sample Plots 



♦ >*€»» 


- 1971 


• a 

J*K. 1971 
6.2* 155.13 

7. SO. 7. *0 

801.70 <101.70 


-osih is res. 

6c 12.19 

b C 7 . J S 

701. 57 1C1. 

.C 

10MM IS fiOCH 
iC 17.77 
SC 7.59 

627.60 570. 

.3 

.-ON I h IS *P9IL 

1 2C 27. C2 

12C 7.71 

ccc.c: 1*. 

.iQliIn IS** 

1 5c 28.42 
ISC 7.79 

625.90 15. 

.0 

nilMh IS JUNE 
led 27.07 
UC 7.92 

698.61 15. 

.0 

MONTH IS JULY 
210 27.15 

210 1.1* 

798.29 IS. 

.0 

-On IH IS 6U0. 
2x0 29.32 

2xc 8.98 
848. X2_ 15. 

-Oft !h IS*" SEP7. 

270 11.66 

272 e.98 

972.3* 15. 

.C 

-ON I n IS OCT. 

sac a.Ji 

JCC 9.63 

999.99 15. 

* .2 
- 0 ft I H IS NOW. 
330 8.73 

33c 1C. 92 

9TX.11 15. 

.0 

nOHlh IS DEC. 
16c 3.59 

3 6 C 11.96 

901.90 171. 


*1971 

165.52 

I.J9 

57 701.57 


159. C2 

7.60 

20 999.32 

.3 
1971 
18C. 29 
1.71 

95 18.95 

. 3 

1971 
1 R3.C8 
7.79 

CO 15.00 

*S«71 

213.15 

1.92 

CO 15.00 

‘5»n 

2C9.79 

8.19 

,C3 15.00 

.0 
1971 
2C2.92 

8.98 

,CG __ 15.00 

. *1971 

201.52 

6.98 

,C0 _ 15. DC 

*1971 

167.98 

9.63 

.CO 15.00 

* ? 9 7 1 
168.76 
1C. 92 

,cn is. uo 

.c 

1971 

132.25 

11.97 

.92 92.91 


7.8C 

801.70 

.0 


7.39 
JOi .57 
.C 


7.oC 

276.6 i 

.0 

7.73 
18. XS 

.0 


7.83 

15.00 

.0 

7.99 
15. CO 

.0 


8.2x 

I5.CC 

.0 


8.69 

1S.CC 

.0 


9.2C 
15. c. 

.0 


9.92 
15. CO 
.C 


10. 7X 

15. QO 

.0 


11.58 

30-78 


7 .40 
8C1.7C 
.0 


7. 29 
7CI.51 
.3 


7.t2 

159.26 

.0 


7.60 

18.95 

.0 


7. S3 
15.00 

.0 


8. 15 

15.CC 

.0 


8.51 

15. OC 

.0 


S.C9 

15. OC 

.c 


9.7* 

11. OC 

.0 


1C. 56 
15. OC 

.0 


11. xe . 
15.00 

.0 


11. ss 

X8.01 


7.80 7.80 - 7.82 7.80 7.80 7-80 

801.10 431.70 431.70 40.70 8al.7n 4C1. 

.0 .0 .0 .0 .3 

7.39 7.JS 7.3S 7. 39 7.29 7.39 

701.57 731.57 131.57 7C1.S7 7Q1.57 ?C1. 

.C .0 .3 .0 .3 

7.6b 7.77 7.1s I.U 4.39 6.X6 

97.70 73.88 76.81 114. 6i 221.80 916. 

■0 .0 .0 — - .0 .0 

7.92 4.19 8.5$ 9.25 IC.52 11.22 

18. X5 18. XS 18.95 18.45 10.91 126. 

.3 .g .0 .0 .3 

8.12 8.50 9.21 lalss 12.25 17.65 

1S.CC 15. OG 15. QO 15.00 15.00 19. 

■C .0 .C .0 ^ .0 

8.97 9.11 10.31 12.50 16.29 22.31 

* IS.CC 15.CC 15.00 15.00 15.00 *17. 

.0 .0 .0 • C .0 

9.03 1C. 01 11.77 19.71 19.34 23.83 

15.00 IS. CO 15.00 15.00 15.C0 32. 

•C .e .0 - .3 .3 

9.77 11.15 13. 4C 16.93 22.14 29.21 

is. co is. co 15.00 is.ro is.aa 44 . 

.to .to .0 .0 .0 


13.73 
15 . toO 

.0 


11.7b 

15.00 

.0 


12.79 

15.00 

.0 


11.96 

116.19 


12.92 15. C7 16.9* 22-62 22.98 

15.CC 15.30 IS. 30 rt.xS x25. 

.3 .0 .0 .0 


13.67 16.95 16. x5 19.71 14. T6 

15. OC 15. 3D 20.81 J-x.77 871. 

.0 .2 .0 .0 

19.23 19.23 19.23 14.23 lx .23 

23.9b 430.29 979.11 97x.ll 97x. 

.0 .0 .3 .0 


7.8C 7.80 

70 801.70 SCI. 70 

.0 .3 


7.39 7.39 

, S 7 701.57 701.57 

.3 .0 

8.5b 6.65 

.75 59*.37 627. bC 

.C .0 


11.9*; 11.66 

19 947.51 tea. ci 

.0 .0 

18.93 18.69 

22 390.56 625.40 

.0 .0 


22.69 23.16 

12 9XC.35 696.61 

.0 .0 


24.17 29.39 

60 622.68 798.29 

.0 .0 

24.96 29.59 

53 771.19 494.92 

.0 .0 

23.07 23.12 

C9 926.5 7 972.39 

.0 .0 


19.77 19.77 

99 996.95 999.94 

.0 .0 


19.23 19.23 

11 979.11 974.11 

.0 .0 


11.96 

931.40’ 


11.9b 11.56 11. 46 11.96 

901.40 9C1.HC 901. *C 931 


11.46 11.96 

• X 0 901. XO 901.40 


Figure 5. Sample output - Lake Keowee, 1971 
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TEMPERATURE PROFILES FOR LRKE KEQWEE ’ 1971. 

(DEPTH IS KERSURZO FROM THE DEEPEST POINT OF THE LRKE) 

1 STRT IONS 50Q-SGS) 
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Fiaure 6. Sample plots - measured average temperature profiles (Stations 
SO 0-506) vs predicted temperature profiles. Lake Keowee, 1971 
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TRR7 1 F I CRT 1 QN CYCLE FOR LRKE KEQMEE 1971-i979 


Solid Lines {No Discharge) 

Broken Lines (Discharge - Mid-layer Temperatures) 


equilibrium Temperatures (a) 

( » Surface Temporaturos(j|t) 

V i\ [\ 


Mid-layer Temperatures (a) 



v, U 




1912 1913 . 1914 1 915 • ISIS \f 19 1 
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Figure 7. Sample plot 
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APPENDIX B 

FORTRAN PROGRAM LISTING 
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KJ 

00 


1 

2 

3 

9 

5 

6 

7 

8 
9 


NASA 

C 

C 

C 


C 

C 


syh created on 12 Aug bo at i9:i7;0;> 

ONE DIMENSIONAL MODEL FOR THE -SEASONAL 


THERHOCLINE 


DIMENSION T 120) ,AV (20 1 ,CBC 20) , ZI201, A I 2DI , XKZ1 20 ) .ROW 1201 
DIMENSION 0MI20I , T2(2Q ) , X TDD ( 10.360) 

DIMENSION DELTEM|12).0Pl 12) 

CHARACTERED MONTHS 112 1 

J=1 ,121 /*JAN.* , *F£&.*, *HARCH*, * APRIL* . * MAY 
SEP T.* , *OCT .*, * NOV • * , * DEC • * / 

STORE RESJLTS ON MA6NEHC TAPE READ JRJN=1 


10 


DATAIHONTHSIJ) ,. 

11 


C* JULY* , * AUC. * , * 

12 

C 


13 

L 


1*» 

C 

IF YOJ NEED TO 

15 

C 

OTHERWISE JRUN- 

16 

C 


17 


READ 1 t JRUN 

la 


RE AO 1 ,IYEAR ,02 

19 


READ 1 ,PI , A1 . A 2 

2n 


RE AD 1 s A5 -A6.T0 

21 


READ 1 ,R6 ,Rt iRB 

22 

1 

FORMAT ( 1 

23 


HMI-U 

29 


1 11 )=0. 

25 


J 1 M-l 

26 


TDD-O. 

27 


0 V£ -D . 

28 


CALL AREAS (A) 

29 


J -1 

3D 


JU = 1 

3] 


J J-0 

32 


NDAYS=0 

33 


N 0 AY S 1 -0 

39 


TIME=0. 

35 


TlMEl=U. 

36 


time?=o. . 

37 


TIME3-0. 

38 


T1KE9 -0. 

39 


TE-TO 

90 


DO 20 1=1 ,12 

81 


T II )=T0 

92 


T2 11 )=T0 

93 

2D 

CONI IN JE 

99 


DO 22 1=2,11 

95 


2 II )=DZ/2. «1 


TNI20) 


.•JUNE*. 


K» 

«£» 


46 

47 
46 

49 

50 
LI 
5? 

53 

54 

55 

56 

57 
56 

59 

60 
61 
62 

63 

64 

65 

66 
67 
66 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 


22 CONTINUE 

Z (1 2 1 =H 

DT=(0.4*0Z**2l/100U.U 
QP2=579.07363*!60.**2l *24 • 

CALL Y EARS (SELTEM, uQPP.lYEARl 
CALL CCU (QP.DELTEH, IYEAR.OTl 
N=0 

0KE6A=2.*PI/36Q. 

T 112>=T0 
T 12 “TO 
J TOT =1 
M J-l 

ROW (12 1=(A4*A5*T(121*A6*(T( 121 1** 21 ♦62. 4 
ROkCP-ROU (121 *CP 

CALL EuUlLl UN ,T£, XK, TDEtt • *IND ,HSGL1 
33 IFIMJ.EQ.l J0 ELTm2=oM( 1 1 ~T ( 71 

FRVt'L = (R6 ♦R7*S1MOHE6A*TIME*R0 11**2 
XKZ0=(R9*R10*SIN(0KEGA*TIME*R8 i 1 
AV(1 1 =A1 ♦ A2* (I(ll-4.1*A3*(T( ll-4.1**2 
XKZ (1 1-XK20* ( l *S1 GHA*A\j! ( 1 1*G* ( IH-21 11 1**21* 

1(3. *7(1 l*T (31^4. *T (211 /(2.*0Z*FRVELI 1**(N-11 
00 SO 1=2 ,11 

AV ( 1 ) - A1 *A2* (T( 11 “4.1 ♦ A3* ( T( I 1-4.1** 2 

XrZ (I 1=XKZ0* (1*SICKA*AV(I 1*G* 1(H-2(1 1 1**21* 

1 (1 (I *1 1-1 (1-1 1 1/ (D2*FRVEL1 1** (N-ll 

ROW (I 1=(A4 «A5*T III ♦ A6* (1(1) 1**21*62.4 
90 CONTINUE 

R Ow (12 I - ( A4 *A5*T ( 1 21* A 6*4 T (121 !*♦ 21*62.4 
AV (12 1 =A1 ♦ A2* (T( 12 1-4 . 1*A3*( T( 121 -4. 1 **2 
XKZ (12 l=XK/0* ll*SI£»MA*AV( 12I*G*: (H-Z( 1 21 1 **21 * 

1 (3. *T (11 1*1 (91-4. *1 1101 1 /( 1 «5*02*F RVEL1 l**(l«-l 1 

ROwCP=ROW (121* CP 

CALL SMOOT h (X KZ, XK2U , X K2L, NDAY SI, TN12, 1 1 2 , T, DT 1 , 0 Z 1 


79 

60 

61 

62 

63 

64 

65 

66 

67 

68 
69 

90 

91 

92 
9 3 

94 

95 


902 DO 969 1=1,12 

IF (XKZ (I 1 .LT.XKZL1 XKZ( Il-XKZL 
IF (XKZ (I 1. CT.XKZ01 xKZl I1=XKZ0 
989 CONTINUE 

DO 91 1=2,11 
F 1 =DT/ (ROW (II * CP* A ( ll 1 

F 2 = 1 (ROU (I l*ROw(I« 11 1/2.* (A (1 1 *A( 1*11 1 /2. 

1*1 XK Z ( I 1 ♦XKZ (1*111 /2.*(T( 1*11 -1(11 l-( (ROul 11 
2* ROU (I — 1 1 1*IA ( II ♦ A (1*1 1 1/4 ( XKZ ( 1 1 

3* XKZ (I — 1 1 1./2 • * (T ( 1 1— T ( I— 1 1111 / (0Z**21 
IF (IYEAR.LL.19731DlLTH 2=0.D 
IF (IYEAR.LE.197310 P2=0 . 0 
777 F3zRGw (I 1* LEL TEM(J»l*CP*UP(JUl 
F 31 =R0 « (I 1*0ELTM2*LP*0P2/A III 

F4i = (ROU ( 1 1*CP*GP2 /( 1 .5*021 1*0 EL TM 2* ( T(I*11-T( 1-1 1 I 

F 4 = ( RO «f (I )*CP*QP(JhJ / ( 1 ,5*0Z1 1 *DEL1EM( JU1*(T1 1*1 l-I (1-111 

IF (T (1*1 l.LE.T ( 1-1 >lF4 = ([R0w(I l*CP*viP( Jkl/( l.S*LZl )*OELTEM (JUl* 



96 


IF II 11 *1 I.LE.T I I-l )IF41=(R0W(I )*CP *QP 2/C 1 . 5*02 1 )*DELTK2 

97 


F5-T II 1 

98 


F6zfl.5* lEXPI-0.75* IH-ZlI 1 1 » )*IHSOL) 

99 


F7--0.75*All) 

ICO. 


Ffiz-F6*F7 

101 


TO zT |8 ) ♦DELTEM I JU ) 

10? 


IF (I .LT.8 IXAKzq. 

103 


IF (I .6E.8 IXAKzl. 

104 


IFITII l.GT.TGlXHrO. 

105 


IF IT ll).L£.TD)XM=i. 

106 


TD2=0ELTH2 ♦TISl 

107 


IF 11 II l.GT .ID2JXM1-0. 

108 


IFIT 11 I.LE.T 02 )X Ml -1. 

1G9 


IFIT.LE.51XTK=1. 

no 


IFlI .6 T.S IXTK-1 • 

111 


TNI I 1 = IF2«F3*XAK*XH*F4*XMl*F41*F31*XTMF8)*Fl*F5 

11? 

91 

CONlINJE 

1 13 


TNI11-TI2J 

1 14 


Th-ITN 112 1 ♦TDEWI/2 .0 

115 


FW=9.2 *0.4 6* (U 1N0* *2 1 

116 


BETA=0.35*0.015*TM*0.00129ITH**2I 

117 


X K - 1 4 . S *0 . 05* T N 1 12 1 ♦BETA^F **0 . 47*F k) * 4 .2 32 *1 9. /5. » 

1 18 


TE-10EW*HS0L/XK 

119 


C0NS1 -1 1 . 5*X K*D? ) /|K OuC P*XKZ 1 12) I 

120 


TE1 1 -TN 111) 

121 


TE10=TNI10 ) 

122 


SHE AT = (ROWCP*DELTEK( JU)*UPI JW ) 1/lAI 12)*XK) 

123 


IFITO.GT .1NI12) )G(j TO 14 

124 


GO TO 15 

125 

14 

TNI 12 1 = I4.*TNI11) -TNI 1 0) ♦CONS 1*TE*5HE A T*CONS 1 1/I3.*C0NS1 ) 

126 


GO TO 16 

127 

15 

TNI 12 )-l4.*TNIll)”TNI 10) ♦CONS 1*TE 1/ 1 3. ♦CONS 1) 

128 

16 

TS-TNI12 ) 

129 


CALL M I X I T ITN .A) 

1 30 


T I ME-T IME *01 

121 


TIME2-TIME2-*DT 

1 32 


T IME3-TIME3*DT 

133 


TlME4zTlHE4*DT 

134 


TlMESzTlMES^Of 

1 35 


DO 929 I - 1 ,12 

1 36 

\ 

T 2 II lzTN'iI) 

137 

929 

CONTINUE 

138 


T 1 2 — T 1121 

1 39 

600 

T N12--T N 1 1 2 ) 

140 

601 

00 92 I-l ,12 

141 


T 11 l-TNII ) 

142 

92 

CONTINUE 

143 


JZJ^I 

144 


TIMEI-TIMEI+DT 

145 


IF INDAYS • GE. 360 ) T IME 3-T1 HE 3-360. 0 



1 fit 
197 
1 *18 
1 *«9 
1 53 

151 

152 

153 
15*» 

155 

156 

157 


IF (NDAYS.LE.36017 IML2=TI ME2-360. U 
IF CNDAYS. CE.3601T i ME - TI ME - 36 0 . 0 
IF CNDAYS. GE. 360 IT iKL*l=Tl ML <1-360.0 
IF <NDAYS.CE. 360 ) T IMF 5=TI ME5 -360 . 0 
IF CNDAYS. GE. 360) JJ=0 
IF < NO AYS. 6E. 360 I J.-l 
IF flYEAR. GT.1979I 60 TO 99 
IF CKO AYS. GE. 360 ) I YEAR — IYLAR^ 1 

IF CNDAYS.GE.360ICALL CCW < uP, DEL Tt H. I YEAR , DTI 
IF (NDAYS . GE. 3601 CALL YiEA RS < SEL7EM ,GCPP , I YEAR 1 
IF CNDAYS. GE. 360) JT0T=JT0T«1 
IF (NDAYS. GE. 360} J1M=JIM* 1 


158 

159 

160 
161 
162 
163 
16*1 

165 

166 

167 

168 

169 

170 

171 

172 

173 
1 7*1 
1 75 

176 

177 

178 

179 
160 
lfcl 
162 
163 
18*1 
165 
1 c6 

167 

168 
169 
190 


1*/1 

192 

193 


501 


502 


66 


98 


988 


313 

700 


I F f TIM E8 • GE. 1 . 0 ) C C TO 501 
GO TO i02 
H MI -MM 1*1 
XTOD|JlH f MMI)=TD 
TlME9=TIME*l-l . 

CONTINUE 

IF <NDAYS.GE.360)N0AYS = 0 

DO 66 1=2.10 

CBII | = (T(l*l|-TflJ)/15. 

CONTINUE 

CB <1 1 = <T (2 l-T 111 1/7.5 
CBCll I = <T <121-1(11 J 1/7 .5 
I F (T IKE1.GE.3U . 1 Go TO 98 
TDD=TDD*TO 
DVt =DVE*1 . 

GO TO 33 

NDAYS=TIME2 

TDG=TUU/DVE 

PRINT 988 ( (CP ( J»’ J) , J. J-l | 12) 
FORMA! ( 1 X , 12F 1 0. 1 ) 

TIME*! =0 • 

M Ml — G 
J J = J J * l 


JW^J-^1 
N0AYS1 -T IM£3 
M J t"M J * 1 


D ELTM2 -OH <M J) — T ( 5 I 
IF IHJ. uE. 12 1H J=1 
CONTINUE 
DO 700 1=1 ,12 
T <1 1=TN<I ) 

IF (JRUN.E0.21 GO TO 111 

CALL STORE <T,AV r CB,2, A,XKZ» ROU.TN.OM, T MONTHS .T2 ,GP . 
£S P, 5 I S'2 A,Ri » Rt4 *«5,k6 f R7,R6 , R9, RIO, wP 2 » FRL UtL , ROW CP , of » 
CXK/iO .TE.NDAYS ,IN12 ,Tl2,Fl,F2,F3,F3l t F*!l,F5 t F6 f F7,F£,TD t TD2 
CNDAYS1 ,TlME,TlME2 t 1 IME 3 ,1 YEAR , MJ , XK , TOG , J 1 


» 



I 


w 

KJ 


194 

195 

196 

197 
193 
199 
2 GO 
2G1 
2G2 
203 
2 04 

205 

206 
207 
2G3 
2 09 
210 
211 
212 

213 

214 

215 

216 

217 

218 


111 CONTINUE 

CALL E0UIL1 (TN , T E , XK, TDEW , WIND ,H SOL 1 
PRINT 920 .MONTHS (JJt, IYEAR 

920 FORMAT (2X t 'MONTH IS • t 2X , Afc t 2X , 14 1 
PRINT 101 .(.DAYS ,TE ,Xn 

101 FORMAT (IX ,16 ,2F9. 2 I 

UR1TEI6.9I N0AYS,ITCll t I=l,121 
URITEI6.71 XK20,(XKZm, 1-1,121 

IF I (IYEAR ,E0. 1973. AND s , NDA Y :> .G E .21 01. OR. (I YEAR. GT. 1973 ) 1 
CURITEI6.18 )TDO (DELIEMI JU-1 ) 

18 FORM AT (IX > ‘THE AVERAGE MONTHLY OISCH. TEMP. = • t F5.2,5X. 
C* OELT A-T = * t F5 • 2 1 
12 FORHAT (IX .11F10.21 

9 FORMAT fix ,16. 12F9. 21 

7 FORMAT (1X.13F9. 21 

TlMEl=TlMEl-30.0 
T 0 0 -0 . 

ove=o. 

IF (IYEAR. GT. 1979 1GG TO 99 
GO TO 33 

99 PRINT 921 ,J 

921 FORMAT (2X. 'TOTAL NJMbER OF COMP JTA TION S = * »1 15 . * X 12*1 
E NO FILE fi 

STOP 

ENO 



1 

2 

2 

4 

5 

6 

7 

8 
7 

10 

11 

12 

13 

14 

15 

16 
17 
13 

19 

20 
,21 

22 


APE A S SYM CREATES ON 12 AJG 80 M 13 : 05:27 
c this SUBROUTINE CONTAINS he areas of 
C A OOP AIN (LAKE KEO.EE), AT TWELVE 

C HORIZONTAL CROSS-SECTIONS. 

C 

C 


SUBROUTINE AREAS(A) 
DIMENSION A (20 1 
AC 0 NS= 1 Q .**8 
A(l) =Q. 032 S*AC 0 NS 
A ( 2 )= 0 . 055 *AC 0 NS 
A ( 3 IS 0 . 20 C BACONS 
A ( 4 ) rO » 55 Q*A CONS 
A ( 5 )= 1 . 125 *AC 0 NS 
A (6 ) =1 • 8 * A CON S 
A ( 7 ) Z2 . 575 »a CONS 
A ( 6 ) - 3 . 55 * ACONS 
A (9 J = 4 .70 ACONS 
A ( 10 )= 5 . 825 *A CONS 
A(ll )= 7 . 25 *AC 0 NS 
A (12 )=fl. 00 d*A CONS 

return 

END 


)h POOR QUALITY 


33 


ecu SYM CREATED ON 12 AU 5 SO AT 13:00:09 


1 

C 

THIS SUBROUTINE CONTAINS THE C CN DENSER 

? 

C 

COOLING WATER. ASSUMES THAT COMPUTATIO 

T 

Ik 

c 

r 

START IN 1971. 

5 

Sm 

c 


6 

c 


7 


SUBROUTINE CCU (QP, uELTEM, I YEA 9,0 T> 
DIMENSION UP ( 1 2 ) ,u£LTE M( 12 1 

3 


9 


IF (IYEAR.GT.1979)GQ TO 1 1 

10 


IYE A-IYEAR-197Q 

11 


A COS T -10 • 0 

15 

1 

GO TO ( 1 * 1 *3*4 ,5,6, 7,8,9) , I YEA 
00 10 1=1 ,12 

14 


OP (I 1=0.0 

15 

in 

DELTEM (I >=C.O 

16 


GO TO 11 

n 

3 

DO 12 1=1,6 

13 


OP II >=0.0 • 

19 

12 

OELTEM (I >=L.O 

20 


OP 17 ) = 1890 »2*AC0ST 

21 


QP (8 ) -19 1 0 . 39 A CO ST 

22' 


OP 19 ) =2 1 70 .7* A CO ST 

23 

* 

OP (10 >=2232.59AC0S T 

24 


OP (1 1 > =217C.7*AC0S T 

25 


OP 112 ) =3284.6»AC0S T 

26 


OELTEM (7 > = 5.3 

27 


OELTEM (3 >=4.b 

2.3 


OELTEM (9 > = 5.3 

29 

30 


DELTEM (10 >=7.3 
DELT EM (11 >=7.7 

31 


OELTEM (12 >=4.1 

32 


GO TO 11 

33 

4 

QP (1 > =3069 • 3* A COST 

34 


OP (2 >=3069 .49AC0ST 

35 


QP (3 > =2976 . 9* A CO S T 

36 


QP (4 >=230 7 .3* A CO ST 

37 


QP (5 >=2164 .6*AC0ST 

33 


OP (6 >=4171 .8*AC0ST 

39 


OP (7 >=S334 .69 ACOST 

40 


QP (8 )=4727.1*AC0ST 

41 


QP 19 > =596 1 . 4* A CO ST 

42 


OP (10 > =49 5 3.49AC0S T 

43 


QP (1 1 > =4202. 1 9 ACOS T 

44 


QP (12 > =5225.6*AC0S T 

45 


DELTEM (1 )=4.2 

46 


DELTEM (2 >=7.4 

47 


OELTEM (3 >=8.4 

48 


DELTEM (4 ) =8.Q 

49 


OELTEM (5 )=2.7 

50 


OELTEM (6 > = 6.0 

51 


DELTEM (7 > = 5.0 

52 


OELTEM (8 > = 4.8 

S3 


OELTEM (9 >=5.8 

54 


OELTEM (10 > = 3.5 

55 


OELTEM (111=7.9 

56 


OELTEM (12 >=5.9 

57 


GO TO 11 

58 

5 

OP (1 >=4612 .49AC0ST 

59 


QP (2 > = 3694 .99AC0ST 

60 


QP (3 >=5456.89AC0ST 

61 


OP (4 )=557Q ,89A COST 

62 


QP (5 >=6494 .39AC0ST 

63 


QP (6 >=6574 .29AC0ST 

64 


OP (7 >=7104 .2*AC0ST 

65 


OP (8 >=7510 .19ACQST 

66 


QP (9 >=7201 .6* A CO ST 

67 


OP (1C >=6993.4*AC0ST 


34 


63 


OP 11 1 ) =7467.1*AC0S T 

69 


OP (121 =685C.9*AC0S T 

70 


DELTEM (1 >=6.2 

71 


DELTEM (2 >=4.o 

72 


OELTEM <3 >-c.2 

73 


OELTEM (4 ) = c • 3 

74 


DELTEM (5 >=6.8 

75 


OELTEM (6 >=6.6 

76 


OELTEM (7 ) =6 • 3 

77 


OELTEM 18 >=7.8 

73 


OELTEM <9 > = 7.4 

79 


OELTEM (10 >=7.7 

60 


OELTEM (11 >=6.5 

81 


OELTEM (12 )=9.4 

82 


GO TO 11 

83 

6 

3 P (1 >=6069 .3*ACQST 

o4 


OP (2 >=4440 .2*AC0ST 

65 


OP (3 >=4874 .3*AC0ST 

36 


OP (4 >=4272 • 1 * A CO S T 

67 


CP (5 >=3970 .7*AC0ST 

63 


OP (6 >=5197 .6 * A CO ST 

69 


OP (7 )=5830 .G*AC0ST 

9 0 


QP (8 1=7248 ,3*AC0ST 

91 


QP (9 )=6785 .4*AC0ST 

92 


OP (10 >=5637.8*ACOS T 

9? 


OP (1 1 ) =58Q9.2*AC0S T 

94 


OP (12 ) =4914.8*AC0S T 

95 


OELTEM <1 > = 10.6 

96 


0 E L T EM ( 2 > = 7. 3 

97 


OELTEM (3 1=7.1 

93 


DELTEM (4 >=5.1 

99 


DELTEM (5 > = 5.8 

ICO 


OELTEM (6 >=9.3 

1C1 


DELTEM (7 > = 7.4 

1 02 


OELTEM (8 )=o.»5 

1C3 


0 ELT EM (9 > = B.Q 

1C4 


DELTEM (10 ) =7 . 8 

105 


OELTEM (1 1 >=6.7 

1C6 


DELTEM (12 )=8.4 

1C7 


GO TO 11 

1 03 

7 

OP (1 >=5045 ,8*AC0ST 

1C9 


OP (2 >=4985 .2* A CO ST 

110 


OP (3 >=S113.S*AC0ST 

111 


OP (4 > =6013 »6*ACQST 

1 12 


QP (5 ) =6 302 .4*A CO ST 

i n 


OP (6 ) =4 3 9 5 • 3*A CO S T 

1 14 


QP (7 > =50 3 8 • 6* A CO S T 

1 15 


QP (8 >=5708 .9*AC0ST 

116 


OP (9 ) =6964 . 0* A CO S T 

117 


OP (10 ) =6754.7*AC0S T 

118 


OP (1 1 > =469 7.6* ACOS T 

119 


OP (12 )=5354.6*ACQS T 

120 


DELTEM (1 > = 12.5 

121 


OELTEM (2 > = 11.4 

122 


DELTEM (3 > = 10.4 

123 


OELTEM (4 > = 11.4 

124 


OELTEM (5 >=9.4 

125 


OELTEM (6 >=8.4 

126 


DELTEM <7 >=7.4 

127 


OELTEM (8 >=5.0 

123 


OELTEM (9 >=5.0 

129 


OELTEM (10 1=3.8 

130 


OELTEM (11 )=6.2 

131 


OELTEM (12 ) =7 • 9 

1 22 


GO TO 11 

133 

a 

OP (1 >=6176 .7* A CO ST 

1 34 


QP (2 >=6444 .6*AC0S7 

135 


QP (3 >=5195 ,7*AC0ST 

1 36 


OP (4 >=4811 ,8*AC0ST 

1 37 


QP (5 >=4984 ,2*AC0ST 

1 33 


OP (6 >=5659 .9* A CO ST 

139 


OP (7 >=7053 .S*ACGST 
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t 


* 


14C 


OP <a )=7914.9*AC0ST 

141 


QP <9 )=6557.3*ACOST 

142 


QP (10 > =7407.4*ACDS T 

143 


QP (11 > =6065.1»AC0ST 

144 


OP <12 ) =65D J.5+AC0S T 

145 


DELTEM <1 > = 9 .C 

146 


DELTEM (21=11.0 

147 


DELTEM (31=13.2 

1 48 


DELTEM (4 )=9.7 

149 


DELTEM (5 ) = 1G. 1 

150 


DELTEM (6 )=b.l 

151 


DELTEM (7 > = 7.9 

152 


DELTEM (fi ) =7 . 5 

153 


DELTEM <9 > = 7.6 

154 


DELTEM (10 >=6.2 

155 


DELT EM (1 1 >=8.4 

156 


DELTEM (12 ) =7. 2 

157 


GO TO 11 

153 

9 

QP (1 >=7207.7*AC0ST 

159 


OP (2 ) =73 1 9 .9*A COST 

160 


QP (3 >=74 19 .5 * A CO ST 

161 


QP (4 >=7275 • 8 » * A C 0 5 T 

1 62 


QP (5 >=418 9 .i*ACOST 

163 


QP (6 >=533 1 .2* A CO ST 

•1 64 


QP (7 >=4733 .3*AC0ST 

1 65 


QP (8 >=4 733 • 3* A CO S T 

166 


QP (9 )=4 733.3*AC0ST 

167 


QP (10 ) =4733.3*AC0ST 

163 


QP (1 1 ) =4733.3*AC05 T 

169 


QP (121=4733. 3*AC0ST 

1 70 


DELT EM (1 > = 10.3 

1 71 


DELTEM <2 1=10.4 

1 72 


DELTEM (3 1=9.6 

1 73 


DELTEM (4 )=9.9 

174 


DELTEM 15 ) =c • 2 * 

175 


OELT EM (6 >=7.1 

176 


DELTEM (7 > =5 • Q 

1 77 


DELTEM (8 >=5.0 

178 


DELTEM (9 ) =5 • 0 

179 


DELTEM (10 >=5.0 

130 


DELTEM (11 )=5.0 

1S1 


DELT EM (12 ) =5.0 

182 

11 

RETURN 

163 


ENO • 
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w 


1 

2 

3 

9 

5 

6 

7 

8 
9 

10 

11 

12 

13 

19 

15 

16 

17 

18 

19 

20 


EQJIL1 SYM CREATED ON 11 JuN 80 AT 11:00:00 

SUBROUTINE EOOILK TN.TE.XK, TDEk.X TN t XTE, XXK,WlND v HSOL 1 
DIMENSION TK ( 2 0 ) ( XTN(201 
2 RE ADI5 ,1 ITOEtf .WIND ,HSOL 

1 FORMAT (1 

U I NO -k I N D* U * 9 £» 

H SOL -HS0L*3.6855 

T M - 1 T N 1 1 2 I ♦TDEUJ /2 .0 

FU-9 .2 ♦O.'Jfa* IU IND**21 

BET A-G • 35 *0.01 5#TM *0.0012* ( TH**2> 

XK— 9.5 -*Q«0 5*TW 1121 ♦BET A+FU + O. 97 *’» k 

XK:XK*<4.232* (9./S. 1 

TE-TDEU*HSCL/XK 

X TM = 4 X TN 1 1 2 1 ♦ T OE W 1 /2.0 

XFk-9.2^0.96* I W IND**2 1 

X BETA=C.35*0.015*XTM*0 . 00 1 2*4 X TM«*21 

XXK=9.5*0.U5*XTN4l 2 1 ♦ XBET A *XF U + 0 ,97*XFk 

XXK-XXK*9 .232* I9./5. J 

XTE-TDEU*HSOL/XXK 

R£TURN 

END 


1 

2 

3 

4 

5 

6 

7 

8 
9 

15 

11 

12 

13 

14 

15 

16 
17 
13 

19 

20 
21 
22 

23 

24 

25 

26 
27 
23 

29 

30 

31 
12 

33 

34 

35 

36 

37 
33 
39 
*40 

41 

42 

43 

44 

45 

46 

47 
4 8 

49 

50 

51 

52 

53 

54 

55 

56 

57 
53 
59 
cO 
61 
62 

63 

64 

65 

66 
67 
63 
69 


INPJT SYM CREATED ON 12 AJG 80 AT 13:01:13 

3.0. 6.6 9,167.0 
C., 9 . 3,264. 4 

6.3. 5.23 .264.4 
7.5, 6.72,457. 

17.2 .7.5.430.5 
18* 3,5.65 ,473. 

iQ . ,o.48,4C9. * 

19.44.5 .75 .428.2 
13.33,5.77,329. 

13. 88,7 .G2 ,261.3 
2. 88, 7*53^24 7.7 

5.5. 6.3.147.7 
1.67 ,6,69 ,173. 

-2.22,9.26 ,257.6 

I. 11 ,9.20,352.5 

6.67 .6. 72 ,448. 

II. 1 1,7 .53 ,433.6 
1 3. 1 3,7 .95 ,564.3 
13. 7 7,6 .64 ,493.3 

22.22.6 .07 ,453.5 

18.3.5.47.336.3 

11. 5 .7. 17 ,298.1 

5. 9. 7.1 3.220.9 
4. ,6. 8, 146. 

1. . 7.22 .162.7 
—1. , 7.3 ,279.5 

10. . 7.1 .348.5 
7. 7, 8.4 4,449.3 

14.3 ,6.S5 ,449.5 

20.25.3 .Q4 ,507.7 

22.2 ,5. 32 ,496.9 

21.7.5.1.391.6 
20. 3 ,6. 3C3 ,333.4 

13.5.7.1.341.7 

7.2.8.14.247.6 
3. 2, 5. 6 ,15 4. 

5.2.5.3 ,191.4 
C., 5. 8, 226. 9 

6.3. 7.7.326.1 
in. 7,8. 73 ,397.7 

17.2,6.3,43b. • 

17. 8 ,6. 98 ,559.3 
21. , 5.2 ,459.5 
21 . , 5.8 7,480. 

1 7. 5 ,6. 74 ,339.2 

10.2 ,5. 7,302.5 

6.0. 7.2 .231.1 

3. 8. 6. 9 .181.9 

3.0. 6.393.191.4 

3.5 . 7.6 14 .226.9 

2.2. 9.6 .326.1 

7.2. 7.6 .397.7 
17.5 ,4.6 ,436. 

19.0 ,5.82 ,559.3 

21.3 ,5. 10,459.5 

21.0. 5.4.480.8 

16.2 ,7. 3,339.3 

12.4 ,7. 7,3G2 .5 

7.9.6.9.231.1 
2. Or 7*2 ,131.9 
-1.0,7.4,209.8 

3.2. 8.5 .310.9 

3.9. 7.9.338.6 

11.2.7.6.496.9 

14.0. 7.3.448.4 

13.3.6.4.480.2 
19. 8 ,5. 9,483 . 3 

13.0 ,6.65 , 4 8 C • 4 

15.4 ,7. 13 ,345.1 
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70 

8.2 , 7.2 1 ,287. 

71 

1.0, 7.27,237. 

72 

-1.5 ,8.2,195. 

73 

— 6.6 ,8.04 ,205 

74 

-2. 75,8 .4 ,317 

7*5 

6.0, 7.7 ,328.5 

76 

10.2 ,7.6,427. 

77 

15.4 ,6.2|473. 

78 

19. ,6.7 ,543.3 

79 

20.2 ,5. 8 ,551 . 

80 

2 0.7, 5.4,423. 

81 

IS. 7,5. 3,350. 

82 

9.2, 7.2,256.6 

83 

7.0, 7.5 ,196.2 

84 

0.4 , 7.2 ,178.2 

85 

-2. P ,7. 9,227 . 

86 

-5.0 ,6.8 ,3C8 . 

to 7 

1.2, 7.6,408. 

6 8 

9.6, 7.6 ,429. 

89 

14. ,6.7 ,513. 

90 

1 9.4 ,4. 7,593 . 

91 

20.8,5.7,568. 

92 

2n. 8 ,5. 1,461. 

V 3 

15.5 ,5.7,365. 

94 

9.3 , 6.6 , 369, 

’95 

9.0 , 5.8 ,232. 

96 

0.4, 7.3 ,191. 

97 

— 3. 33,8 . o , 2 Q o 

93 

C.G , 7.2 ,25 1 . 

99 

b.O , 7.9 , 373. 

ICO 

9.2 , 7.6 ,479. 

1 u 1 

14. ,0.7 ,513. 

1C2 

1 9. 4 ,4.7,598 . 

1G3 

20. 8 ,5.7,568 . 

104 

20. 3 ,5. 1 ,461 . 

1C5 

1 5.5 ,5. 7,387 . 

1 C6 

9.3,6.6,369. 

1C7 

9.0 , 5.8,232. 

1C3 

u. 4, 7.3,191. 

109 

0.4 *7.3,191. 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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4 

5 

6 

7 

8 
9 

10 

11 

12 

II 

14 

15 

16 

17 

18 

19 

20 
21 
22 


MIXIT SYM CREATEC ON 12 AUG 80 AT 13:26:57 

C 

C 

C THIS SUBROUTINE MIXES STABILIZES UNSTABLE 

C TEMPERATURE PROFILES. 

C 

c 

SUBROUTINE MIX IT(TN,A) 

• DIMENSION TN ( 20 ) ,A IZO) 

ICO 00 10 1=1 ,11 

I F ITS ( I *1 ).GE.TN( I»»GO TO 1 

IF t (TN II )-TN(I*l I ) .LT .0.0)60 TO 1 

T A V = ( T N ( I ♦ l ) + T N ( I) )/2. 

TN (I ♦! >=T A\i 
TN (I ) = TAV 
1 CONTINUE 

10 CONTINUE 

TMAX =AM AX 1 (TN ( 1 > ,T M2 T , TN< 3 > , TN(4) , TN ( 5) 

C , T N ( 6 ) ,TN ( 7 ) , T N { 3 ) , T N 
C ( 9 ) ,TN (10 ) ,TN(ll),TN(12)) 

IF (TN ( 12 ) .LT. TMAX) 00 TO lOu 
300 RETURN 

END 
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1 

2 

3 

4 

5 

6 

7 

8 
9 

IQ 

11 

1? 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 
21 
3? 

33 

34 

35 


38 

39 

40 

41 
4? 

43 

44 

45 


PLOTTER SYM CREATED ON 12 A J6 80 AT 12:56:46 
PARAMETER Nri4 ,NN- 12 , NTI ME - 12 , ND-1J 0 
DIMENSION IBuFIlOOCl 

DIMENSION TEMP c SO I ,DEEP(501 , TEtf PSlNO 1 i DEEPS (ND I (OP (NN ) 
DIMENSION T IN) ,AV( M ,CB(N1 ,Z(N),XKZ(N) ,TEG INDI ,THF CN01 
< ? W * N * f T M N 1 .DMINI, T2CN1 . A ( N I , ZED (NO ) 
DIMENSION LI 15 01 ,E^(50),E3l501,E4IS0),E6(50l,E5(50;, 
CE7 (50) ,ED(50) 

CHARACTER*6 MONTHS IN) 

CHARACTER*o IBCD 
M-l 
L to 


C 

C 

C 

C 

C 

c 

100 


READ JRUN=1 IF YOU DESIRE PLOTS FOk MEASURED DATA 
READ JHUN-2 IF YOU DO MOT. 

NOTE : IF PLOTS FOk SEVEN STATIONS ARE NOT 

AVAILABLE, LINES 35 TO 46 MUST faE MODIFIED 
READ 100.JRUN , JYEAR 
FORMAT ( ) 

I COUNT tq 
X Z0=G. 

JO TQ 

CALL PLOTS (IBUF, 1000, 11) 

CALL PLOT (0.0 ,7.0, -3) 

DO 1 I -1 ,N T IH £ 

^ V ,C o , 2 , A ,|X KZ , R 0 J , TN , DM , TZ, MONTHS, T2 ,QP , 
CCP .SIGMA ,R3 ,R4 ,R5,R6 t R7,R8 ,R9, R1 j, iP2, FREVLL.ROUCP.OT . 
CXK20,TL,NDAYS,f»Jl2,Tl2,Fl,f2,f 3,Fh,FKl,F5,F6 # F7 t F& # T6. 
CNDAYS1 JTIML.TIME2, TlME3,IYLAR,MJ,XK,TDD,,;i 
I COUNT =f COUNT *1 

if (I Count. lt . 96 i go to 3 33 

IF (JRJN.EO .2) GO TO 200 

RE AD (5 ,8 1 (DEEP ( INK 1, TE MPl I NK ) , INK- 1 , NS TOP ) 

DO 15 KL -1 ,50 

,8 1 uEEPIKL) ,E 1 ( KL) , L 2( K L 1 , E 3( KL ) , E 4 (KL 1 ,E5 (KL I , 

Ct u (KL 1 |l 7 I K l_ 1 

n^rn »6)AEl,6El , C E 1 , D E 1 , E L 1 , F E 1 , L E 1 , HE 1 , 0 E 1 
DLLP 4KL 1-AE1 * 

El (KL 1 rfiEl 
E 2 (KL ) TCP 1 
E 3 (KL l-DEl 
E4 (KL)tEE 1 
E5 (KL ) TFE1 
EG (KL )=GE1 
E7 (KL l thEI 


,TZ (NN) 
,TSU (NO) 


TD2 , 



•p 

K» 


46 


47 


48 


49 


50 


51 


5? 

15 

53 

200 

54 

16 

55 


56 


57 


58 


59 


60 

222 

61 

223 

62 

201 

63 


64 


65 


t6 


67 

9 

68 


69 


70 


71 


72 

202 

73 

a 

74 

333 

75 


76 


77 


78 


79 


80 


ei 


62 


63 


64 


65 


66 


67 


68 


89 


90 


91 


9? 


93 

203 

94 

444 

9 5 



ED IKL l=OEl 

IFIE3IKL1.EQ. 0.0160 TO 16 
IF IOEEPIKL l.tO. C-l .1 1G0 TO 16 

TEMP IKL 1 = 1 El IKLI *E 21 KL > «E3 I KL > «E 4 1 KL 1 ♦ L5 1 K L 1 ^E bIKL I « 
CE7 IKL ) 1/ED IKL 1 
TEMP IKL1=E3IKL1 
CONTINUE 
CONTINUE 
NST0P=KL-1 

IF IJRUN.EQ.2IG0 TO 201 

00 222 JIJ=1,5Q 

IF (DEEPIKL l.EO.I-1 .1 100 TO 223 

RE ADIS f 8 I AEl.BEl ,C El ,tJEl , EE 1, FE1 , GE1 , HE1 , 0E1 

IF IAE1 .EQ, l-i I 160 10 223 

CONTINUE 

CONTINUE 

CONTINUE 

CONS 2 — 1 • /O • 304 8 

IFtJRUN.E0.2lG0 TO 202 

00 9 INK -1 ,NS TOP 

OECP UNK )=E0NS2*DEEPl INK! 

DEEP I INK I =150. -DEEP I INK! 

DEEP INSTOP ♦11=0.0 

DEEP INSTOP *2 1 = 2 INN 1/1 .5 

TEMP INSTOP ♦11=0.0 

TEMP INSTOP + 2 ) =30.0 /l.S 

CONTINUE 

FORMAT I 1 

JG=JOM 

L =L^1 

T S J |L 1 =T 112 1 

X 2D=XZ 0*30 . 

ZEDIL 1 =XZO 
TEMPS IL 1 = TEMP 111 

1 COIL 1=TE 

THE IL l.= IT 17 1 ♦ T I 6 1 ) /2. 

I BCD=MONTHS I JO 1 

z Inn ♦ 1 1=0.0 
Z INN+2 1=2 INN1/1 .5 
T INN +1 1=0.0 
T INN *2 1=30. /I. 5 

CALL AXIS 10.0 ,0.0, 6HTEMP.IC1 , -8,l.t»,0.0,Tll3l,Tll411 

CALL AXIS 10.0 ,0.0, 9HDEPTHIFT1 , 9, 1 . S, 90 .0,7 113 I , 2(14 1 1 

CALL FLlNElT, / ,-NN , 1,0,01 • ' 

IFIIC0UNT.GT.961G0 TO 444 

IE IJRUN.E0.21 GO TO 203 

CALL DASHL ITEMP ,DEEP, NSTOP, 11 

CONTINUE 

CALL SYMBOL! 1 .0,0.5,0.14 , 1BCD, 0.0, 61 
CALL PLOT 12.25 ,0.0, -31 



4 = 

U 


96 

97 

98 

99 
ICQ 
101 
10 ? 
103 
109 

105 

106 
1C7 
108 
109 

11 ? 

112 

113 

114 

115 
1 16 
1 17 
118 
1 19 
120 
121 
122 

123 

124 

125 

126 

127 

128 
129 
1 20 

131 

132 
1 33 
134 
1 35 
1 36 

137 

138 

139 


IF I JO. EQ.4 .OR, J0.EU.8 1G0 TO 3 
GO TO 1 

3 CALL PLOT 1-9.0 ,-2. 25,-31 

1 CONTINUE 

CALL PLOT 1-2.25,0.0,-31 

CALL SYMBOL 1-6.75, t». 75 , .14 ,4 1H TEMPERA T LRE PROFILES FOR LAKE KEOWEE 
C ,0.0,91) 

P1=JY£AR 
M Y - J Y E AR 

CALL NUMBER 19 9 9 . ,9 99 . , 0 . 14 , PI , 0. U . 01 

CALL SYMB0L1-6.75, 6.5 , 0 . 1 , S4H I DEPTH IS MEASURED FROM THE DEEPEST P 
COINT OF THE L A KE 1 . 0.0 .54 l 
CALL PLOT 18.0 ,-9.25,-31 
PRINT 2, MY 

2 FORMAT tlX . 'THE PLOTS FOR*, 15,* ARE COMPLETE*) 

IF IM. £0.9 ICO TO 6 

M = M*1 

J YEAR rJY EAR ♦! 

GO TO 5 

6 CALL PLOT (6.0 ,0.0, -3)' 

DO 13 1=1,96 
13 DEEPS(I)=ZEO(I) 

DEEPS (97 1=0.0 
DEEPS (98 1=3240.0/9.0 
TSU (1091=0.0 
TSU(llQ)=35./5. 

TEO ( 1 09 1=0.0 
TEQ(llCl)=35./5. 

T HF ( 109 I =0 .0 
THF (1101=35. /5. 

TEMPS (97 1=0.0 
TEMPS (98 1=35. /5. 

ZE 0 ( 109 1=0.0 
ZED(UQ)=3240./9. 

CALL PLOT (0.0 ,2.0, -3) 

CALL AXIS (0.0 ,0.0, 12HT IME 1 N DAYS, -12, 9.0, 0.0, ZED (109 ) ,ZED (1 10) 1 
^C ALL i AXIS (0.0 ,0.0, 16HTEMPE RA TURES (Cl , 16, 5.0, 9u. , TSU (109 1 ,TSU 

CALL FLlNE (7EU ,TSJ, -10 6,1.2.11 1 
CALL FLINE (ZED, TEO ,108, 1,2, 2) 

CALL FLINE (ZED, THF, -108, 1,2,01 
IF (JRUN.EQ.21G0 TO 204 
CALL UASHL (DEEPS, TEMPS, 96, 11 
204 CONTINUE 



1*10 
HI 
1 «*2 
1 *43 
194 
1*15 
1*16 
1*)7 
1*48 
1*)9 


152 

153 

1 5*» 


C 

C 

C 

C CHANGE TITLES TO SUIT NEEDS < *4 LINES1 

C 

C 

CALL SYMBOL CO.O, 6.0,0. 1*4. 

C46HSTR ATIF1CA1 ION CYCLE FOR LAKE KEOkEE I 9 71 -1 979 . 0.0 ,96 1 
CALL SYMBOL (0*0 ,0. 10*0 .10* B7H 1971 1972 1973 

C 1975 1976 1977 1978 1979.0.0*87) 

URITEC6.7I 

7 format fix. *all plots are ngu complete* ,//, • normal job exit*) 

CALL PLOT 115.0 ,0.0 ,-3) 

STOP 

f ti r. 


-c 

-fc 




1979 


-fc 

Ul 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

II 

14 

15 

16 

17 

18 

19 

20 
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READER SYM 
C 
C 
C 

c 
c 
c 


CREATED GN 12 AUG 80 AT 13:21:45 


THIS SUBROUTINE READS THE MAGNETIC TAPE 
CONTAINING THE COMPUTED RESULTS. 


SUBROUTINE READER I T ,AV ,CB . 2 . A , XK2 , ROM, Tn, OM.TZ, MONTHS ,T 2.QP. 
CCP,S1GMA,R3,R4 ,RS.R6,R7,R8.ft9,RlO. CP2.FREVEL.R0lfCP.Dr, 

CXK20 .TE.NDAYS .TN12 , T 1 2 . F 1 , F 2. F 3,F 3 1 , F 4 1 , F 5 , Ffc . F7 . Ffi . T 0 ,TD2 , 
CNDAYSl *T I M E *T IME2 , TIME* » I YEAR, MJ , XK, TDO.J. NCASE, SF.EDEPT ,V0L1 
DIMENSION T (201 ,AV (20 1 »CB ( 20) *2(201 , A ( 20 1 , XK 2 ( 20 1 1 
CROb (20 l.TN 120) ,DMi201 , TZ,’2Q1, T2f 201.QPU21 
CHAR ACTER*6 MONT HS (121 
CONTINUE 

READ (6 t END=l) CTl I J1 , I J= 1 , 121 ,1 A V(U 1 ,1 J=1 .12 I . 

CJCBdJ I, I J = 1.121 . (ZdJl.lJ-l, 121, (AdJl.IJrl.I2l. 

C ( XKZ II J1 ,IJ=1 ,12), IROW(IJ) . U=l, 121, I TN(IJ1.1J=1. 121, 

C(DH(U 1 ,1 J-1 ,12) ,( TZdJl.ljri, 121', (MONTHS! I J1 , Ijrl.12 1 , 

C ( T2 ( I J I ,lJri ,121 . 

CIQP IIJ 1 ,XJri ,121 , 

CCP,SIGHA,R3,R4jR 5, R6.R7.R8 , R9 , RIO, CP2.FRE VEL » ROW CP , DT , 
CXK20,TE,NDAYS,lM2,Tl2.Fl,F2,F3,F31,F41,F5,F6.F7 i F8,TO,T02j 
CNDAYSl .TIME ,T IME2, TIME 3 , I YEAR ,MJ , XK, TOD, J .NCASE ,SF .EOEPT , VOL 

return 

ENO 
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28 

29 
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SMOOTH SYH CREATED ON 12 AUG 80 AT 14:34:30 
C 

C THIS SUBROUTINE CORRECTS THE EDDY UIFFDSIVITY 

C IF VARIABLE T IHE STEP IS REQUIRED, • D Tl • SHOULD 

C BE CHANGED TO *DT * INI THE CALLING PROGRAM. 

C 

SUBROUTINE SHOOTH4XKZ,XKZU p XKZL,NDAYSl t TN12,T12.T,OTl ,D2l 
DIMENSION XK2 420) , 1 4201 
DO 93 I =1 , 12 

IF 4XK2 41 ). CT.XKZUI XK24I»=XK2U 
IF 4XK2 41 l.Ll.XKZL) XKZ4I)=XKZL 
93 CONTINUE 
NEL=0 

DO 96 1-2 • 12 
IF 4XK2 (I ) . EG.XK2L) AEM-I 
96 CONTINUE 

IF 4NEW.EQ.OIGO 10*77 
DO 55 1=1 , NEU 
XKZ Cl )=XKZL 
55 CONTINUE 
77 CONTINUE 

IF4NDAYSl.LE.6O.OR.NDAYSl.6T.30OI GO TO 29 
IF 4TN12.GE.T12 I GO 10 19 
IF 4TN12.LT. T12 ) GO TO 39 

19 X MIN =AMlN 1 4XKZ 4 1 ) , XKZ 421,XK24 3),XKZ44) , XK 2 4 5 1 , XK Z 4 6 I , X KZ 4 7 ) , X KZ 

146 1.XK249 ) ,XKZ410> V XKZ<111 ,XK?(12> > 

DO 82 1=1 ,12 

IF (XKZ4I1.EQ.XMIN1 GO TO 61 
82 CONTINUE 
GO TO 29 
81 I H I N =1 

DO 70 1=1 , 1M1N 
XKZ 41 I =XKZ 4 1M IN ) 

70 CONT IN JE 

GO TO 29 

39 XMAX=AMAX1 4XKZ 4 1 1 , XKZ 4 2 ) , XKZ4 3>,XKZ44) , XKZ 4 51 , XKZ 46 1 ,XKZ ( 7 1 ,X K2 
1 46 ) ,XKZ 49 ) ,XKZ 4 10) ,XKZ ( 1 1 ) V XKZ 4 1 2) ) 

DO 62 1 = 1 .12 

IF 4XKZ 4 I l.EO.XHA*) ISO TO 61 
62 CONTINUE 
GO TO 29 
61 I H AX =1 

DO 5G I =1 M AX , 1 2 
XKZ4I1=XKZ 4IMAX) 

50 CONTINjE 
29 CONTINUE 

200 XM AX = All AX 1 4XKZ 4 1 I , X KZ 4 2 I , X KZ 4 31.XK2I91 , XKZ 4 5) , XKZ 46 ) ,XKZ 47 ) ,XKZ 

146 1 ,XKZ 49 ) ,XKZ 4 10) ,XK2 Cl 1 ) f XKZ 4 12) ) 

011 = 441.9*02**2 )/XMAX 

RE IJRN 

END 
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STORE 

C 

C 

C 

c 

c 

c 


SYH CREATED ON 12 AJG 60 AT 13:19:47 

THIS SUBROUTINE STORES THE COMPUTED RESULTS ON 
HA6NETIC TAPE. 


™°£r UTINi: , ST0PECT ■*V,CB,2 t A,XK2 t ROH, TN.DM »T2, MONTHS »T 2 »QP , 
CCP t SlbMA v R3 ( R4 t R5 , R6 , R7 , R6 , R9 , R 1 0 . OP 2 , FRE VEL , R OW CP , OT , 

CXK20 ,T E.NDAYS , TN 1 2 ,T1 2 .F 1 . F 2. F 3. F 31, F 4 1 , F 5, F6. F7 . F£ , TO ,TD2 , 

,ME2 » f JHEitlYEAfe.MJ, XK.TOD.J.NCaIe.If.^OEPT t V^Ll 
DIMENSION T (20). AV (211) ,C6 ( 201 ,2 ( !:0) , A ( 20 i * XK 2( 20 1 * * 

R0WI20 I.TNI201 .DH120), 12(20), T2I2G), 

a ri • « . w w 


CR 
COP 112 1 


CHAR ACTFR*6 MONTHS (12) 



CTQPIIJ ) ,ljrl,12| , 

£S£l£ I §£* , £ 3 »S‘‘7 R ?A H6 ! R7 * R8 A R9 J R1U » CP 2 , FRE VtL . RGUCP . DT . 

CXKZO .TE.NOAYS t IN 1 2 , T 1 2 ,F 1 . t 2, t 3 , F Si , F 1 1 , F 5 ,fI, FT, F & ,1 & ,7 D2 , 

C ?°i Y ri. , I I o E,T lME2 » 1 IHE3.IYEAR.MJ. XK, TOD, J, NCASt ,SF, EOEPT , VOL 
t fiy I ILl O 
RETURN 
END 


onnn 


YEARS SYM CREATED ON 12 AUG SO AT 1J:1U:C3 

THIS SUBROUTINE PRINTS THE YEAR TITLE. 


99 


SUBROUTINE Y E A R S < S E LTEt M, G 0 PP, I YEA » ) 
PRINT 99.IYEAR 

FORH AT t59X ,17 ( •* M , /, S9X, *♦*, 1 SX , •+• 
C* *• ,2X , ‘YEAR = • ,I‘+,2X,»*» , /*59X, 
C./.59X ,17 ( '*») ) 

RETURN 

END 


/ » 5 9 X , 

1 5X, * *♦ 


